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In this paper we revisit dynamics of a quantum open system. First of all we explain our 
purpose in a short manner. See as a general introduction to this subject. We consider 
a quantum open system S coupled to the environment E. Then the total system S + E is 
described by the Hamiltonian 

Hs+E = Hs®1e + ^s®He + Hj 

where Hs-, He are respectively the Hamiltonians of the system and environment, and Hi is the 
Hamiltonian of the interaction. 

Then under several assumptions (see pQ) the reduced dynamics of the system (which is not 
unitary !) is given by the Master Equation 

^^p = -r[Hs,p\-V{p) (1) 

with the dissipator being the usual Lindblad form 

V{p) = \Y. {^]^^P + P^]^^- - 2^.P^]) • (2) 

Here p = p{t) is the density operator (or matrix) of the system. 

Similarly, the equation of quantum damped harmonic oscillator (see [T], Section 3.4.6) is 
given by 

—p = —i[ua^a, p] — ^ {a)ap + pa) a — 2apa^) — — (aa^p + paa^ — 2a^pa) , (3) 

where a and are the annihilation and creation operators of the system (for example, an 
electro-magnetic field mode in a cavity), and p, v are some constants depending on the system 
(for example, a damping rate of the cavity mode). 

Since this is one of fundamental equations in quantum theory it is very important to con- 
struct the general solution. In [1] or [2] some methods to construct a solution are presented. 
However, a direct and clear method deriving a general solution has not been given as far as we 
know. It may be difficult to treat because ([H]) is an operator equation. In fact, in [2]; section 
18.3.1 it is stated that "it is extremely difficult to solve the master equation directly". However, 
we can solve the equation (I3l) completely in the operator algebra level. 
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First we separate the equation ([3]) into two parts and consider the one coming from the 
Lindblad form 

—p = — ^ (a^ap + pa^a — 2apa^) — — (aa^ p + paa) — 2a) pa^ . (4) 

Interesting enough, we can solve this equation completely by use of some method in [3], [Ij. 
Let us rewrite (jll) to be more convenient by use of the number operator = a^a ([a, a^] = 1) 

(5) 



d 



+ + + '^Z AT AT \ P — 1^ 

—p = papa^ + ua^pa ^i^P + pN + p) + P 



where we have used aa^ = + 1. Note that [A^, a] = —a, [N, a''] = 

In order to solve the equation we can use the method in [5] once more. For that we review 
a matrix representation of a and on the usual Fock space 



^ = Vectc{|0),|l),|2),|3),---}; \n) 



like 



1 
V2 

V3 




|0) 





1 
V2 





= e-'%^ (6) 



N = a^a = b^b 



V 



where e*^ is some phase. Note that of course [b, b''] = 1. 
For a matrix X = (xij) G M(jF) 

/ 

Xii Xi2 Xi3 
X21 X22 X23 



X 



X31 X32 X33 



(7) 



we correspond to the vector X E j^dimc:^' 

X = (xij) — > X = (xii, Xi2, Xis, ■ ■ ■ ; X21, X22, X23, ■ ■ ■ ; 3:31, X32, X33, ■ ■ ■ ; (8) 

where T means the transpose. The following formula 

AXB = {A®B^)X (9) 

holds for A,B,X e M{T). 
Then ([5]) becomes 



— p = <j pa ® {a^Y + z/a^ ® a - ^-^(^ ®l + l®iV+l®l) + ^-^1 ® l| P 

^ih^h + uh^ - ^^-^{N ®l + l®iV + l(8)l) + ^^-^1 ® l| P, (10) 
so that the solution is formally given by 

By the way, from the old lesson in [3] we know a method to calculate ( iTTl) explicitly. By 
setting 

i^3 = ^(iV®l + l®A^ + l®l), K+ = b^®b\ K_ = h®h, (12) 
we can show the relations 

[i^3, = K+, [K-s, K^] = -K^, K^] = -2/^3 

easily. Namely we have the sm(1, 1) algebra. The equation (fTTl) can be written simply as 

p{t) = e^V{-^++M^-(M+-)i^3}^(0), (13) 
so we have only to calculate the term 

Now the disentangling formula in [3] is helpful in calculating fll4l) . We assume that there 
exists a unitary representation (group homomorphism) p : SU{1,1) C SL{2;C) — > U{T-C) 
where Ti. is some Hilbert space related to the Fock space JF above, and 

dp{k+) = dp{k_) = K., dp{h) = 



where {k+, k_, ks} are the generators of 1) algebra 



oi\ /oo\ l/lO\ , 

\, k^=\ \, k, = -\ ]■ k^^kl (15) 

00/ \-10/ ^\0-l/ 

It is very easy to check the relarions 

[ks, k+] = k+, [ks, k-] = -k-, [k+, k_] = -2/c3. 

From (dH) 

^t{uK++fiK--{fi+u)Ks} ^ ^t{udp{k+)+fi.dp{k-)-{fj.+u)dp{k3)} ^ ^dp{t{vk++fik--{p+v)kz)) 

where 



e*^ = exp { t 



co^h (^t) - sinh (^t) ^ sinh (^t) 

si^h (i^t) cosh (^t) + ^ sinh (^t) 

The Gauss decomposition formula 



c d] I 1 i I J U 1 



; ad — be = 1 



gives 



tA 



f 
f 



-sh(i^i) + ^sinh(ii^t) 





cosh sinh (^t) ) 



X 
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Moreover 



tA 

e = exp 





cosh I ^-^ 





\ 



X 



exp 



yo ; 

^ - log (cosh (^t) + g^sinh(^t)) 



V 
/ 







exp 



\ 



log (cosh (^t) + ^sinh(^t)) ^ 



X 



^ sinhl 







V cosh(4^I) + ii±^sinh(ii^ 

^ sinh (^t) 

^sinh(^t) 



/ ^ sinh 

= '"P|^cosh(4^t)+^±;^s 

-2 log 



■A; . X 



exp 



cosh , 

^ 2 



( ^^inh(^t) 
'"P|^cosh(^t)+;^sinh(^t)'- 



X 



by ffT^ . Since p is a 



group homomorphism [p[XY) = p{X)p(Y)) and the formula p{e^^^ = 
= k+,k3,k_ (see ( lT6i) ) we have 



/ ^ sinh (^t) 



X 



exp 



-2 log (cosh [f^t\ + ^i±^smh (^iZLil()) <i,(t3)' 



^sinh(^t) 

sinh (^t) 



X 



-dpikJ) 



^^Pl^cosh(4^t) + 
As a result we have the disentangling formula 

/ JiL. sinh (^f") 

' P 1^ cosh (^t) + ^ Sinh (^t) ^ 

/ ^sinhf^ 



X 



exp ( —2 log ( cosh ( ^^—^ 



p-u^ 



lE. sinh (^t) 



X 



'^P|,cosh(4^t) + ^sinh(i^t) 



(17) 
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In the following we set for simplicity 

^ sinh (^t) ^ sinh (l^t) 



(^t) + Sinh (^t) ' ^ ^ cosh (^t) + ^ Sinh (^^t) 
.(,^cosh(/i^.).;i±^s.nh(^.). (IS 
Therefore ffTTj) becomes 

p(t) = e"^* exp (^(t)^^" ® 6^) exp (- log{F{t)){N ® 1 + 1 ® iV + 1 ® 1)) exp ® 6) p(0) 

under (|T2l) . By use of a = e*^6 and = e~*^6^ in ([6]), and by some calculation 



e 2 



p(t) =^exp (G(t)at®a^) jexp (- log(F(t))iV) ® exp (- log(F(t))iV)' 



X 



exp {E{t)a ® (a^)^) p(0) (19) 
where we have used A^^ = A^. Coming back to matrix form by use of we obtain 

^ ^ n=0 

:^a-p(0)(at)-l exp(-log(F(t))iV)]a^ (20) 



ml 

^ m=0 



This is indeed complicated ! 



Next we treat the full equation ^ 

—p = —iuj{a)ap — pa^a) — ^ (a^ap + pa)a — 2apa^) — — (aa) p + paa) — 2a^ pa) . 
This can be rewritten as 

^ p = I -iiuKo + i^K+ + pK_ - {p + i^)K3 + ^^-^1 (g) l\ p (21) 



dt 

in terms of Kq = N^1 — l^N. Then it is easy to see 

[Ko, K+] = [Ko, Ks] = [Ko, K_] = (22) 

from f|T2l) . That is, operators arising from the Hamiltonian of harmonic oscillator and its 
Lindblad form commute. Therefore 

^^^^ ^ ^-iiotKo^t{i^K++tiK--{ti+u)K3 + l^im}-^(^Q^ 

= e'^'exp {-iujtKo) exp {G{t)K+) exp (-2 \og{F{t))Ks) exp {E{t)K^) p(0) 
= e'^'exp {G{t)K+) exp {{-lootKo - 2\og{F{t))K3}) exp {E{t)K_) p(0), 



so the general solution of the equation that we are looking for is just given by 



n=0 

J2 ^^^a^p{0){a^r \ exp {{tcot - log(F(t))}iV)}a" (23) 



. 771=0 



by use of (fTQl). 

If z/ = 0, we have 



^ m=0 



I ^'1 —^ty^'' I 

p(t) = e-(^+-)*^ \ aXO)(at)"- }> e"(^-*-)*^. (24) 



To construct the general solution to the equation Q is very important in not only Physics 
but also Mathematics, and we have finished it. 

In this paper we revisited the quantum damped harmonic oscillator and constructed the 
general solution in the operator algebra level. 

The model is very important to understand several phenomena related to quantum open 
systems, so the general solution is needed. We finished it. 

On the other hand we are studying some related topics from a different point of view, see 
[6] and [7j. However, we make no comment on them in the paper and will report some results 
in another paper. 

Lastly, we conclude the paper by stating our motivation. We are studying a quantum 
computation (computer) based on Cavity QED (see [S] and [^), so in order to construct a more 
realistic model of (robust) quantum computer we have to study severe problems coming from 
decoherence. See also [10], [11] for example. This is our future task. 



Appendix 



In this appendix we explain our method in terms of the classical damped harmonic oscillator, 
see for example [6] and [7|. 
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The differential equation is given by 

X + 270; + u^x = (7 > 0) 



(25) 



where x = x{t), x = dx/dt and the mass is set to 1 for simplicity. In the following we treat 
only the case a; > 7 (the case a; = 7 may be interesting). By setting y = x the equation is 
rewritten as 

^ I nr \ I n 1 \ I nr \ 

(26) 



d 


( x\ 




~dt ' 


\ y 1 


H 





1 ^ 


M 


f x\ 


w2 -27 I 




[ y 1 





f -7 1 








[ -1 ) 




\ -7 / 



Noting 

^ 1 

-uj^ -27 

(see ( IT5l) ) the general solution is given by 

x(0) 
Z/(0) 



-7I2 + k+ + LU^k^ + 27^3 



xit) 

y{t) 
( 



g-7tgt(fc+WA:_+27fc3) 



-It 



Q.O'i,{\/ ijP' — '-fH) + 



\ 



7 sin(-^ cij^— 7^t) 



0)^—7^ 



cos( Y^cu^ — 7^^) 



7 sin{-Y/i^^— 7^t) 



cj^ —7^ 



x(0) 

y(o) 



.(27) 



For the details of calculation see [6]. 
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